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Abstract 

The correlation dimension and limit capacity serve theoretically as lower and upper bounds, respectively, of 
the fractal dimension of attractors of dynamic systems. In this paper, we show that estimates of the correlation 
dimension grow rapidly with increasing noise level in the time-series, while estimates of the limit capacity 
remain relatively unaffected. It is therefore proposed that the limit capacity be used in studies of noisy data, 
despite its heavier computational requirements. 

An analysis of Singapore wind data with the limit capacity estimate revealed a surprisingly low dimension 
(-2.5). It is suggested that further studies be made with comprehensive equatorial weather data. 

PACS number: 47.52.+j, 47.53.+n, 92.60.Wc 

1. Introduction 

Estimations of the fractal dimension of meteorological time-series serve a fundamental role in 
the development of dynamic models of meteorological phenomena. These dimensional estimates 
provide bounds for the number of independent variables necessary to model the system, and assists 
one in determining the appropriateness of a model. 

There have been numerous efforts to determine the dimension of weather and climatic attractors 
[1,2,3], however the equatorial weather system has yet to be subjected to a thorough analysis. 
Previous dimensional estimates have been based upon the use of the correlation dimension, but the 
reliability of this measure is questionable as it is unstable when perturbed by noise in the observational 
time series. 

In this study, the stability of the estimates of the limit capacity and the correlation dimension 
is investigated by considering the Lorenz attractor perturbed by various levels of noise. 

The limit capacity dimension of a 22-year equatorial wind time-series, with a sampling 
frequency of once a day is then estimated. 

2. Basic Concepts 

The phase portrait of a dynamical system can be reconstructed from the observation of a single 
variable by the method of delays as proposed by Takens [4]. Takens proved that almost all smooth 
<i-dimensional manifolds can be embedded in a (m=2fi?+7J-dimensional space while preserving 
geometrical invariants. The observational time seriesx^,;^,...^^ represented by the set of vectors 
x={x t jc t+i ,... yC t+(m . 1)T } (where t-l,2,...,N-[m-l]T) in the reconstructed phase space. 
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The delay X is usually chosen for which the autocorrelation function, 

C^) = —Y^(x-x)(x i+% -x) (1) 

(where x is the arithmetic mean) is close to zero, thus minimizing the statistical dependence among 
the coordinates of the vectors. 

In practice, one does not know a priori the dimension of the dynamical system, and the 
embedding dimension which is necessary for the phase space reconstruction. Thus, the dimensional 
estimate is computed for increasing embedding dimensions until the dimensional estimate stabilizes. 



2.1. Correlation Dimension 



The correlation dimension theoretically provides a lower bound to the fractal dimension of an 
attractor, thus it satisfies the following inequality 

d c < d f (2) 

The correlation dimension is obtained by considering the cumulative correlation function, 
defined by Grassberger and Procaccia [5] as 

1 N i n\ 
C(r) = lim — Te r - be, - x . 

N ^~N\% y 11 " W (3) 

where 0is the Heaviside function, such that 9(x) = if x< and G(x) = 1 if x> 0. The Euclidean norm 
utilized in (3) is defined as 



Ir* ~ x j\\ = ^jE( x *+«-i), _ x wh> ) 2 (4) 

The cumulative correlation function is related to the correlation dimension by the power law 

CO)-/' (r-»0) (5) 

This results in the correlation dimension <f being defined as 
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2.2. Limit Capacity 

As opposed to the correlation dimension, the limit capacity, d l serves as an upper bound to the 
fractal dimension of an attractor. Considering this property, and the inequality (2), we obtain 

d c < d f < d, (7) 

Let A^(e) be the minimum number of spheres of radius £ necessary to cover all points of the attractor. 
Takens [4] defined the limit capacity d l as follows: 

A^(e)~e"' / ' (e^O) (8) 

From the relation (8), we obtain the limit capacity d l as 

, ln[jV(e)] 
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2.3. Methods of estimating dimensions 



The reconstruction of an attractor by the method of delays as described in Section 2 is performed 
on the observational time-series. The embedding dimension m which is used in the attractor 
reconstruction begins from the value 2 and proceeds to successively higher dimensions until the 
dimensional estimate stabilizes. For the correlation dimension, C(r) vs r is plotted on a log-log graph, 
and the gradient of the region of the graph which exhibits scaling behaviour yields an estimate of the 
correlation dimension. Similarly, the limit capacity is obtained by plotting N(e) vs 1/e on a log-log 
graph, with the gradient providing an estimate of the limit capacity. 



We consider the effect of noise on the correlation dimension and limit capacity estimates by 
computing the dimensional estimates of the discretely sampled time-series obtained from the x 
component of the Lorenz Attractor, which is described by the following set of ordinary differential 
equations, 



with the following parameters, cr=10, r=28, £=8/3. The Lorenz Attractor is shown in Fig. 1. 

The Lorenz equations were numerically integrated with a fourth-order Runge-Kutta method 
[6], with At=0.005 and x=0.6 (which was the time lag at which the autocorrelation function for the 
discrete time-series of the Lorenz attractor vanished). A discrete time-series with 2048 data points 
of the x component was then generated and subjected to dimensional estimation. Each data point of 
the time-series was perturbed by the introduction of noise with a maximum magnitude which is 1 % 
and 10% of the maximum magnitude of the x component. The noise was generated by a pseudo- 
random multiplicative linear congruential generator proposed by Park and Miller in [7] . 

As can be seen from Fig. 2, the dimensional estimates of the Lorenz attractor without noise 
converge to a finite value with 2.06 for the correlation dimension and 2.31 for the limit capacity at 
an embedding dimension m-(2n+l ), thus satisfying inequality (2.2. 1). The attractor with 1% noise 
as shown in Fig. 3, also behaves similarly with a correlation dimension of 2. 1 83 and a limit capacity 
of 2.342. However, when the time series is perturbed by 10% noise as shown in Fig. 4, inequality (7) 
is not satisfied, as the correlation dimension of 4.678 is greater than the limit capacity of 2.332. 



3. Effect of Noise on Dimensional Estimates 



dx 




dt 



(10) 




Fig. 1 : Lorenz Attractor 



Fig. 2: Dimensional Estimate of Lorenz Attractor 





With 1% Noise 

— ■ Gaussian Noise 

— • - Correlation Dimension 
- * Limit Capacity 














* 





1 2 3 4 5 6 
Embedding Dimension 





With 10% Noise 

— * — Gaussian Noise 

* Correlation Dimension 
- *• Limit Capacity 



















Fig. 3: Dimensional Estimate of Lorenz Attractor 
with 1% Noise 
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Fig. 4: Dimensional Estimate of Lorenz Attractor 
with 10% Noise 



3.1. Analysis of Dimensional Estimates 



The correlation dimension is extremely sensitive to noise, as illustrated by the rate of divergence 
of the correlation dimension from the dimensional estimate of the noiseless attractor upon the 
addition of noise. This is significant, since the increase in the correlation dimension due to noise 
would potentially render the lower bound to the fractal dimension provided by the correlation 
dimension invalid, as in the case of the Lorenz attractor with 10% noise. Thus it would seem that the 
correlation dimension would tend to measure the dimension of the noise as opposed to the underlying 
dynamics which are of interest. A possible explanation for this is that for a noisy attractor, the form 
of the power law relation (5) is 

C (r)~r m r<r 

/?? V / noise 

C (r) ~ /' r>r (11) 

'-'mV / ' noise v/ 

Given that the scaling region for the correlation dimension occurs in the region ( r Q), it follows 
that the dimensional estimate would tend to be higher than that inherent in the dynamics of the system. 

The analogous situation does not occur in the case of the limit capacity as shown from the small 
deviation from the dimensional estimate obtained for the noiseless attractor. An intuitive explanation 
is as follows. In the region of small r, the limit capacity estimate suffers from the effect of saturation 
due to high embedding dimensions and a finite time-series, and hence it is generally necessary to 
obtain the limit capacity estimate at a moderate r region, where the effect of noise is less dominant. 



4. Dimensional Estimate of Wind Attractor 



The time-series analysed are the x and y components of the daily wind velocity observations 
over Singapore. Missing data were linearly interpolated. A Fourier transform was performed upon 
each of the two time-series. 

As shown in Fig. 5 and Fig. 6, there is a conspicuous peak in the yearly cycle range. The long 
term cycles are of such high magnitude that they have to be filtered off to allow for the analysis of 
the dynamics of the weather for periods of less than 140 days. The autocorrelation function is plotted 
for this filtered data and is shown in Fig. 7 and Fig. 8. The value of x for which the autocorrelation 
function vanishes is 4 days for the X-Winds time-series and 8 days for the Y-Winds time-series. 
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Fig. 7: Autocorrelation of X-Winds 
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Fig. 8: Autocorrelation of Y-Winds 
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Fig. 9: Limit Capacity Estimate of Wind Data 



The limit capacity dimensional estimate was then performed upon the two time-series for 
embedding dimensions from 2 to 10. The dimensional estimate was obtained by plotting the points 
(l/e,N( £)) on a log-log graph. The region satisfying the power law relation would be represented by 
a straight line on the log-log graph. A least-squares fit is then performed to provide the limit capacity 
dimensional estimate. As seen from Fig. 9, the estimated dimension deviates rapidly from the linear 
behaviour which would have been expected if the series consisted of Gaussian noise, thus implying 
that the dynamics are of a deterministic and finite-dimensional nature. It is however essential to 
realize that the choice of range for the scaling region from which the limit capacity dimensional 
estimate was derived is somewhat subjective, and a different range could result in a slightly differing 
dimensional estimate. 



Embedding Dimension 


Gaussian Noise 


X-Winds 


Y-Winds 


2 


2 


1.537 


1.392 


3 


3 


1.939 


1.634 


4 


4 


2.190 


1.901 


5 


5 


2.271 


2.017 


6 


6 


2.350 


2.025 


7 


7 


2.343 


2.109 


8 


8 


2.402 


2.193 


9 


9 


2.430 


2.294 


10 


10 


2.508 


2.313 



Table 1 : Limit Capacity Dimensional Estimate of Wind Data time-series 

As seen from Table 1 , the limit capacity converges to a fractional dimension between two and 
three, 2.508 for the X-Winds time-series, and 2.313 for the Y-Winds time-series. This fractional 
dimension is indicative of the sensitive dependence on initial conditions which characterizes the 
dynamics of the weather system. 

The correlation dimension estimates were also computed, and they exhibited behaviour similar 
to the dimensional estimates for the Lorenz attractor with 10% noise, and did not satisfy inequality 
(7). 

5. Conclusion 

The dimensional estimates obtained from the correlation dimension should be treated with 
caution as the dimensional estimates are not robust when the time-series is perturbed by noise. The 
limit capacity in comparison is less affected by noise and is thus recommended for dimensional 
estimates of noisy dynamical systems. 

The low fractal dimensionality obtained for the equatorial weather attractor suggests that the 
inherent dynamics could be modelled by a low-dimensional deterministic system. The dimensional 
estimate for the Singapore wind data is lower than that reported for temperate systems. This is 
surprising as the equatorial weather systems are usually less organized than the temperate weather 
systems. A comprehensive study with regional data on a range of weather parameters may provide 
a deeper insight into the dynamics. 
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